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We present a general derivation of semi-fermionic representation for generators of 5*17 {N) group 
as a bilinear combination of Fermi operators. The constraints are fulfilled by means of imaginary 

(N : 

Lagrange multipliers. The important case of SU{2) group is discussed. We demonstrate how the 
idea of semi-fermionic representation might be extended to the groups possessing dynamic symme- 
tries. As an example, 50(4) group is considered. We illustrate the application of semi-fermionic 
representations for various problems of strongly correlated physics. 
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INTRODUCTION 



^1 It is known that spin operators satisfy neither Fermi nor Bose commutation relations. For example, the Pauli 
^ ' matrices for S* 1/2 operator commute on different sites and anticommute on the same site. The commutation 
c/2 [ relations for spins are determined by SU{2) algebra, leading to the absence of a Wick theorem for the generators. 
To avoid this difficulty and construct a diagrammatic technique and path integral representation for spin systems 
various approaches have been used. The first class of approaches is based on representation of spins as bilinear 
combination of Fermi or Bose operators [1]- [6], whereas the representations belonging to the second class deal with 
more complex objects like, e.g. the Hubbard [7] and supersymmetric [8] operators, the nonlinear sigma model [9] 
^ , etc. However, in all cases the fundamental problem which is at the heart of the difficulty is the local constraint 
Q ■ problem. To illustrate it, let's consider e.g., first class of representations. Introducing the auxiliary Fermi or Bose 
O i fields makes the dimensionality of the Hilbcrt space, where these operators act, greater than the dimensionality of 
the Hilbert space for the spin operators. As a result, the spurious unphysical states should be excluded from the 
consideration which leads in turn to some restrictions (constraints) on bilinear combinations of Fermi/Bose operators, 
^ , resulting in substantial complication of corresponding rules of the diagrammatic technique. The representations from 
' the second class suffer from the same kind of problem, transformed either into a high nonlinearity of resulting model 
(non-linear sigma model) or hierarchical structure of perturbation series in the absence of Wick theorem (Hubbard 
operators). The exclusion of double occupied and empty states for a S' = 1/2 impurity interacting with conduction 
' electron bath (single impurity Kondo model), is controlled by fictitious chemical potential (Lagrange multiplier) of 
Abrikosov pseudofermions [4]. At the end of calculations this "chemical potential" A should be put A ^ — cx) to 
"freeze out" all unphysical states. In other words, there exists an additional U{1) gauge field which freezes the charge 
-j—; fluctuations associated with this representation. The method works for dilute systems where all the spins can be 
considered independently. Unfortunately, attempts to generalize this technique to the lattice of spins results in the 
^ , replacement of the local constraint (the number of particles on each site is fixed) by the so-called global constraint 
I ' where the number of particles is fixed only on an average for the whole crystal. There is no reason to believe that such 
■ O an approximation is a good starting point for the description of the strongly correlated systems. Another possibility 
^ ■ to treat the local constraint rigorously is based on Majorana fermion representation. In this case fermions are "real" 
and corresponding gauge symmetry is Z2. The difficulty with this representation is mostly related to the physical 
regularization of the fluctuations associated with the discrete symmetry group. 

An alternative approach for spin Hamiltonians, free from local constraint problem, has been proposed in the 
pioneering paper of Popov and Fedotov [10]. Based on the exact fermionic representation for S — 1/2 and 5=1 
5J] ' operators, where the constraint is controlled by purely imaginary Lagrange multipliers, these authors demonstrated 
C3 ' the power and simplification of the corresponding Matsubara diagram technique. The semi-fermionic representation 
(we discuss the meaning of this definition in the course of our paper) used by Popov and Fedotov is neither fermionic, 
nor bosonic, but reflects the fundamental Pauli nature of spins. The goal of this paper is to give a brief introduction to 
a semi-fermionic (SF) approach. A reader can find many useful technical details, discussion of mathematical aspects 
of semi-fermionic representation and its application to various problems in the original papers [10]- [21]. However, we 
reproduce the key steps of important derivations contained in [18], [19] in order to make the reader's job easier. 

The manuscript is organized as follows: in Section I, the general concept of semi-fermions is introduced. We 
begin with the construction of the SF formalism for the fully antisymmetric representation of SU (N) group and the 
fully symmetric SF representation of SU{2) group using the imaginary-time (Matsubara) representation. We show a 
"bridge" between different representations using the simplest example of 5 = 1 in SU (2) and discuss the SF approach 
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for 50(4) group. Finally, we show how to work with semi-fermions in real-time formalism and construct the Schwinger- 
Keldysh technique for SF. In this section, we will mostly follow original papers by the author [11], [18]. The reader 
acquainted with semi-fermionic technique can easily skip this section. In Section II, we illustrate the applications 
of SF formalism for various problems of condensed matter physics, such as ferromagnetic (FM), antiferromagnetic 
(AFM) and resonance valence bond (RVB) instabilities in the Heisenberg model, competition between local and non- 
local correlations in Kondo lattices in the vicinity of magnetic and spin glass critical points and the Kondo effect in 
quantum dots. In the Epilogue, we discuss some open questions and perspectives. 



I. SEMI-FERMIONIC REPRESENTATION 



To begin with, we briefly reproduce the arguments contained in the original paper of Popov and Fcdotov. Let's 
assume first S = 1/2. We denote as H„ the Hamiltonian of spin system. The standard Pauli matrices can be 
represented as bilinear combination of Fermi operators as follows: 

^^.^a]a,-b]h,, <-2a]6„ aj^2h]a,. (1) 

on each site i of the lattice. The partition function of the spin problem Z„ is given by 

Z„ = Tr exp{-l3H„) = i^Tr exp(-/3(^ir + mNp/i'^.p)) (2) 

where Hp is the operator obtained from Ha by the replacement (1) and 

N 

N = Y,{a]ai+b\h^) (3) 

i=i 

(A'^ is the number of sites in the system and /3 = 1/T is inverse temperature). To prove equation (2) we note that the 
trace over the nonphysical states of the z-th site vanishes 

Tr^nphvs eM-P{HF + nrNp/m) = + = (4) 

Thus, the identity (2) holds. The constraint of fixed number of fcrmions Nj = 1, is achieved by means of the purely 
imaginary Lagrange multipliers n = — z7r/(2/3) playing the role of imaginary chemical potentials of fermions. As a 
result, the Green's function 

G = {iojF - e)-i (5) 

is expressed in terms of Matsubara frequencies cop = 2iTT{n + 1/4) corresponding neither Fermi nor Bose statistics. 
For 5 = 1 we adopt the representation of H„ in terms of the 3-component Fermi field: 

crj -^a]a- b]b, a+ ^2 (a] cj + c] bj ) , aj ^ V2(c] aj + 6] Cj ) . (6) 
The partition function is given by 

Z, = Tr{-(3H,) = i-^j Trexp{-P{HF + inNp/m). (7) 

It is easy to note that the states with occupation numbers and 3 cancel each other, whereas states with occupation 
1 and 2 are equivalent due to the particle-hole symmetry and thus can be taken into account on an equal footing by 
proper normalization of the partition function. As a result, the Green's function in the imaginary time representation 
is expressed in terms of up = 2TrT{n + 1/3) frequencies. 

In this section, we show how semi-fermionic (Popov-Fedotov) representation can be derived using the mapping of 
partition function of the spin problem onto the corresponding partition function of the fermionic problem. The cases 
of arbitrary N (even) for SU(N) groups and arbitrary S for SU(2) group are discussed. 
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A. SU(N) group 



We begin with the derivation of SF representation for SU(N) group. The SU(A'') algebra is determined by the 
generators obeying the following commutation relations: 

KiSr]=^iMS!i-€S'.i), (8) 

where a,/3 = 1,...,N. We adopt the definition of the Cartan algebra [22] of the SU(A'') group {Ha} = similar 
to the one used in [23], noting that the diagonal generators S" are not traceless. To ensure a vanishing trace, the 
diagonal generators should only appear in combinations 

JV N 

with ^.s„=0, (9) 

a=l a=l 

which effectively reduce the number of independent diagonal generators to N — 1 and the total number of SU(A'^) 
generators to N'^ — 1. 

In this paper we discuss the representations of SU(N) group determined by rectangular Young Tableau (YT) (see 
[23] and [18] for details) and mostly concentrate on two important cases of the fully asymmetric (one column) YT 
and the fully symmetric (one row) YT. 

The generator may be written as biquadratic form in terms of the Fermi-operators 

5|=^aty^ (10) 

7 

where the "color" index 7 = 1, ric and the ndric + l)/2 constraints 

N 

E««7l«"'''=^^> (11) 



restrict the Hilbert space to the states with m*nc particles and ensure the characteristic symmetry in the color index 
a. Here m corresponds to the number of rows in rectangular Young Tableau whereas ric stands for the number of 
columns. The antisymmetric behavior with respect to a is a direct consequence of the fcrmionic representation. 
Let us consider the partition function for the Hamiltonian, expressed in terms of SU(A^) generators 

Zs = Tr exp(-/3i7s) = Tr' exp(-/3i?F) (12) 

where Tr' denotes the trace taken with constraints (11). As it is shown in [18], the partition function of SU{N) model 
is related to partition function of corresponding fermion model through the following equation: 

Zs = j l[di^{j)P{i^{j))Trexp{-PiHp - /x(j>f)) = / l[dfi{j)P{f,ij))ZF{i^{j)) (13) 

i i 

here P{fij) is a distribution function of imaginary Lagrange multipliers. We calculate P{l^j) explicitely using con- 
straints (11). 

We use the path integral representation of the partition function 

Zs/Zl = y n dix{j)P{ix{j)) exp(^)/ j n dix{j)P{ix{j)) exp(^o) (14) 
3 3 

where the actions A and Ao are determined by 

A = Ao~ / dTHpir), A = V V / drukU, T){dr + l^{j))ak{j, r) (15) 

and the fermionic representation of S\J{N) generators (10) is applied. 

Let us first consider the case ric = 1. We denote the corresponding distribution by P;v^m(/i(j)), where m is 
the number of particles in the SU(A^) orbital, or in other words, 1 < m < N labels the different fundamental 
representations of SU(A''). 
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N 



= X] '^k{j)akU) = m, (16) 



fe=i 



To satisfy this requirement, the minimal set of chemical potentials and the corresponding form of Pjv,m(At(i)) are to 
be derived. 

To derive the distribution function, we use the following identity for the constraint (16) expressed in terms of 
Grassmann variables 



-sm(4n,-m))/sm( ^^^^ '-] (17) 



Substituting this identity into (12) and comparing with (14) one gets 

Pn,M3)) = Iv E«^P - 1) 5{„{3) - n,), (18) 

fe=i ^ ^ 

Mfe = - — (2fc-l). (19) 

Since the Hamiltonian is symmetric under the exchange of particles and holes when the sign of the Lagrange multiplier 
is also changed simultaneously, we can simplify (18) to 



where 



fW,m(M(i)) = sin Urn— — j S{n{j) - /Xfe) (20) 

fc=i ^ ^ 

where [N/2\ denotes the integer part of N/2. As shown below, this is the minimal representation of the distribution 
function corresponding to the minimal set of the discrete imaginary Lagrange multipliers. Another distributions 
function different from (20) can be constructed when the sum is taken from k = N/2 +1 to N. Nevertheless, this DF 
is different from (20) only by the sign of imaginary Lagrange multipliers jl^ = fx^ = — /Zfc and thus is supplementary 
to (20). 

Particidarly interesting for even N is the case when the SU{N) orbital is half-filled, m = N/2. Then all Lagrange 
multipliers carry equal weight 

2, "^/^ 

Pn,n/2{hU)) = ^ E(-l)'^''^ (^(■?') - Z^'^) • (21) 
fe=i 

Taking the limit N oo one may replace the summation in expression (21) in a suitable way by integration. Note, 
that while taking N ^ oo and m ^ oo limits, we nevertheless keep the ratio m/N = 1/2 fixed. Then, the following 
limiting distribution function can be obtained: 

Pn^nMi^U)) ^ exp (-/3M(j)y) (22) 

resulting in the usual continuous representation of the local constraint for the simplest case Tie = 1 

Zs = rr(exp {-(}Hf) S (nj - y )) (23) 

We note the obvious similarity of the limiting DF (22) with the Gibbs canonical distribution provided that the Wick 
rotation from the imaginary axis of the Lagrange multipliers fi to the real axis of energies E is performed and thus 
n{j)N/2 has a meaning of energy. 

Up to now, the representation we discussed was purely fermionic and expressed in terms of usual Grassmann 
variables when the path integral formalism is applied. The only difference from slave fermionic approach is that 
imaginary Lagrange multipliers are introduced to fulfill the constraint. Nevertheless, by making the replacement 

/■ \\ /■ \ { iiTT 2k — \\ _ , . ■. _ N / ii^T 2k — 1\ 
a-k{3,T)) ak{3,T)exp i— — I , ak{j,T) afe(j,r)exp I — — \ (24) 
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we arrive at the generalized Grassmann (semi-fermionic) boundary conditions 



akij,f3) = afc(j, 0)exp in 



2k- 1 

N 



a-kUiP) = afe(j,0)exp -in 



2k -1 
N 



This leads to a temperature diagram technique for the Green's functions 

G"^ij,T)=-{Tra^ij,T)apij,0)) 



(25) 



(26) 



of semi-fermions with Matsubara frequencies different from both Fermi and Bose representations (see Fig. 2). 

The exclusion principle for this case is illustrated on Fig.l, where the S* = 1/2 representation for the first two 
groups SU(2) and SU(4) are shown. The first point to observe is that the spin Hamiltonian does not distinguish 
the n particle and the n hole (or N — n particle) subspace. Eq. (19) shows that the two phase factors exp(/3/in) 
and exp(/3/i(A^ — n)) accompanying these subspaces in Eq. (20) add up to a purely imaginary value within the same 
Lagrange multiplier, and the empty and the fully occupied states are always canceled. In the case of > 4, where 
we have multiple Lagrange multipliers, the distribution function P(/i) linearly combines these imaginary prefactors 
to select out the desired physical subspace with particle number n = m. 



SU(2) 



n=2 




n=4 



n=0 



Re e 




SU(4) 



Re e 



im e-^^^l 

n=l ^-^"""^ 




n=3 


n=4 / 

^ 




\ n=0 










n=2 





FIG. 1. Graphical representation of exclusion principle for SU(Af) semi-fermionic representation with even N , Uc = \ (we 
use p. = inT/2 for SU(2) and = inT/'i, pa = 3i7rr/4 for SU(4)). 

In Fig.l, we note that on each picture, the empty and fully occupied states are canceled in their own unit circle. 
For SU(2) there is a unique chemical potential = ±inT/2 which results in the survival of single occupied states. 
For SU(4) there are two chemical potentials (see also Fig. 2). The cancellation of single and triple occupied states 
is achieved with the help of proper weights for these states in the distribution function whereas the states with the 
occupation number 2 are doubled according to the expression (21). In general, for SU(-/V) group with ric = 1 there 
exists N/2 circles providing the realization of the exclusion principle. 



B. SU(2) group 

We consider now the generalization of the SU(2) algebra for the case of spin S. Here, the most convenient ferniionic 
representation is constructed with the help of a 25' + 1 component Fermi field akU) provided that the generators of 
SU(2) satisfy the following equations: 

s-i S 

5+ = E VS{S+l)-k{k + l)al^,{j)akU), VS{S+l)-k{k-l)al^,{j)ak{j), 

k=-S k=-S+l 

S'^ fc4(j>fc(j) (27) 

k=-S 

such that dirnHp = 2^^^+^ whereas the constraint reads as follows 

k=S 

^3 = E 4(j>fc(j) = ^ = 1 (28) 

k=-S 
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Following the same routine as for SU(iV) generators and using the occupancy condition to have I = 1 (or 25") states 
of the (25"+ 1) states filled, one gets the following distribution function, after using the particle-hole symmetry of the 
Hamiltonian Hs'. 

2i / 2fc-l\ 

P2S+i,i(m(j)) = E sin (^TT^^j - A^fc) (29) 

where the Lagrange multipliers arc /i^ = —iTiT{2k — 1)/(2S' + 1) and k ~ \, [S + 1/2J , similarly to Eq.(19). 




"■+^ "+^(1-^) '^ + ^(1-2^1) ^ 



n n+1/2 n n+1/2 

FIG. 2. The minimal set of Matsubara frequencies for a) SU{N) representation with even N/ SU{2) representation for 
half-integer value of the spin, b) SU{2) representation for integer values of the spin and I = 1. 

In the particular case of the SU(2) model for some chosen values of spin S the distribution functions are given by 
the following expressions 



P2,M))^^S[t^ij) + — ] (30) 



for S = 1/2 



pMj)) = pMj)) = ^ 5 [p^ij) + ^) (31) 



for S = l. 

This result corresponds to the original Popov- Fedotov description restricted to the 5 = 1/2 and S = 1 cases. 




FIG. 3. Graphical representation of exclusion principle for SU(2) semi-fermionic representation for 5* = 1 and 5* = 2. For any 
arbitrary integer value of spin there exists S circle diagrams corresponding to the S different chemical potentials and providing 
the realization of the exclusion principle. 



A limiting distribution function corresponding to Eq. (22) for the constraint condition with arbitrary / is found to 

be 

Poo,M)) ^^eM~mj))- (32) 
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For the case I = m, = N/2 oo and S = {N — l)/2 ^ oo the expression for the Hmiting DF Poo.i{p{j)) coincides 
with (23). We note that in 5 ^ oo (or N ^ oo) hmit, the continuum "chemical potentials" play the role of additional 
U(l) fluctuating fleld whereas for finite S and N they are characterized by fixed and discrete values. 

When S assumes integer values, the minimal fundamental set of Matsubara frequencies is given by the table in 
Fig.2. 

The exclusion principle for SU(2) in the large spin limit can be also understood with the help of Fig.l and Fig. 3. 
One can sec that the empty and the fully occupied states are canceled in each given circle similarly to even-N S\J{N) 
algebra. The particle-hole (PH) symmetry of the representation results in an equivalence of single occupied and 2S 
occupied states whereas all the other states are canceled due to proper weights in the distribution hmction (29). In 
accordance with PH symmetry being preserved for each value of the chemical potential all circle diagrams (see Fig. 3, 
Fig. 5) are invariant with respect to simultaneous change 

^ — M cllld f^particle ^ '^holes' 

C. Prom SU(2) to SO(4) 

We have shown that the general rectangular Young Tableau of size Uc * m is represented hy N * ric component 
fermionic field with Uc diagonal constraints and nc{nc — l)/2 off-diagonal constraints. However, the fully symmetric 
representation (one row) requires only ric + l = 25+1 component field. The general scheme of projected representation 
for SU(N) group is given in [18]. We illustrate this idea on a simple example of 5* = 1. 

We start with 2 * nc = 4 - field representation 

(an, 012, a2i, 022) (33) 
There are two diagonal and two off-diagonal constraints which read as follows: 

ail«ll + «2l'*21 = 1, «I2'^12 + «22<^22 = 1- (34) 

anai2 + 021*22 = 0, ajjaii + 022*21 = (35) 
and generators of SU{2) group are given by 

S~ = Sl = o{i021 -I- 0^2022, = Sf = O21O11 -I- 0I2O12 

25^ = S2 — Sl = 021021 -I- O22O22 — oj^oii — o|20i2 (36) 
Combining definition (36) with constraint (35) we reach the following equations: 

5" = o|i(o2i -h 012) -h (012 + 021)022, S"*" = (oil -h a|2)oii -h 022(012 -h 021), S'^ = 0^2022 - ojiOii (37) 

Therefore, we conclude that the antisymmetric (singlet) combination 012 — 021 docs not enter the expression for spin 
S = 1 operators. Thus, three (out of four) component Fermi-field is sufficient for the description of S = 1 SU (2) 
representation. Defining new fields as follows 

On = /-I, 022 = /i, ■^^'^i^ + 021) = /o, "^("12 - 021) = s. (38) 

where fermions /i, /o, f-i stand for £'■^ = 1,0—1 projections of the triplet state and fermion s determines the singlet 
state, we come to standard S = 1 SU{2) representation (c.f 6) 

5+ = ^/2(/t/-l + /lVo), 5- = \/2(/ljo + 4/i), S, = flfi-fUf-i, (39) 
with the constraint 

rii +no + n^i + Us = 2 (40) 

where = flfa- 

Nevertheless, the constraint (40) transforms to a standard SU(2) 5=1 constraint in both cases = and = 1 
since there is no singlet/triplet mixing allowed by SU(2) algebra. 
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To demonstrate the transformation of the local constraint let's first consider the case Ug = 0. The constraint reads 
as follows 



ni +no + ri-i = 2S 



(41) 



On the other hand, the the states with 25 occupation are equivalent to the states with single occupation due to 
particle-hole symmetry. Thus, the constraint (40) might be written as 



rii + fio + fi-i = 1 

where = 1 — riQ,. The latter case corresponds to Us = 1. 

We start now with definition of 5*0(4) group obeying the following commutation relations 

[Sj,Sk] = iejkiSi, [Pj,Pk] = iejkiSi, [Pj,Sk] = iejkiPi 
where 6 generators of SO (4) group, namely vectors S and P are represented by the matrices 



10 

-^"^ ' 




S- =y/2 



/O \ 
10 
10 

Vo 0/ 



/ 1 

0-10 

\0 



(42) 



(43) 



P+ = V2 



1 


0-10 



With the Casimir operator 



V2 



/O 

-1 

V 1 



S • P = 0, + P2 = 3. 



pz 




-1 

0-100 



(44) 



Unlike SU(2) group, the singlet /triplet transitions are allowed in S0(4) group and determined by P operators. Using 
the definition of singlet /triplet fermions one comes to following representation 



^+ = \/2(/t/_i + /t/o), 5- = ^/2(/lJo + /oVl), = fifi - fUf-i, 



(45) 



P+ = V2{fls - st/_i), P- = V2{s^f, - fUs), = -(4s + st/o). 
with the only constraint 

ni + no + ti-i + ris = 1 
whereas the orthogonality condition is fulfilled automatically. 



(46) 



D. Real-time formalism 

We discuss finally the real-time formalism based on the semi-fermionic representation of SU(iV) generators. This 
approach is necessary for treating the systems out of equilibrium, especially for many component systems describing 
Fermi (Bose) quasiparticles interacting with spins. The real time formalism [24] , [25] provides an alternative approach 
for the analytical continuation method for equilibrium problems allowing direct calculations of correlators whose 
analytical properties as function of many complex arguments can be quite cumbersome. 

To derive the real-time formalism for SU(A'') generators we use the path integral representation along the closed 
time Keldysh contour (see Fig. 4). 
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00 ^ ^ ) + oo 



c. 



t-i p 



FIG. 4. The Koldysh contour going from —oo oo ^ — oo in real time. The boundary conditions on the imaginary time 
segment determine the generalized distribution functions for quasiparticles. 

Following the standard route [26], we can express the partition function of the problem containing SU(A'^) generators 
as a path integral over Grassmann variables tpi = (a(,i(j), ai^M{j))^ where I = 1,2 stands for upper and lower parts 
of the Keldysh contour, respectively, 

Z/Zq = j D4>Djpexp{iA)/ j D4>Djpexp{iAo) (47) 

where the actions A and Ao are taken as an integral along the closed-time contour Ct + Cr which is shown in Fig. 4. 
The contour is closed at t = — oo + ir since exp(— /?i7o) = TV exp ^— J^^ H^dr^ . Wo denote the '0 fields on upper and 
lower sides of the contour Ct as Vi and V2 respectively. The fields ^ stand for the contour Ct- These fields provide 
the matching conditions for ^1^2 and are excluded from the final expressions. Taking into account the semi-fermionic 
boundary conditions for generalized Grassmann fields (25) one gets the matching conditions for '^1^2 aX t = ±00, 



^l,alfe(-Oo) = exp i^^'-^Yf-) V'2,aU(-00), 



V'5',„U(+oo) = V2%ife(+oo) (48) 

for fc = 1, \_N/2\ and a — 1, ...,7V. The correlation functions can be represented as functional derivatives of the 
generating functional 

Z[r]\ =^0^ J D^jDi; exp (^A + ij dt{f]a''^ + ^cr^r?)^ (49) 

where 77 represents sources and the matrix stands for "causal" and "anti-causal" ordcrings along the contour. 

The on-site Green's functions (GF) which are matrices of size 2N x 2N with respect to both Keldysh (lower) and 
spin-color (upper) indices are given by 

To distinguish between imaginary-time (26) and real-time (50) GF's, we use different notations for Green's functions 
in these representations. 

After a standard shift-transformation [26] of the fields ip the Keldysh GF of free semi-fermions assumes the form 

where the retarded and advanced GF's are 

G(«'^'«(e) = (e±i^)-\ /, = /W=)(e), (51) 
with equilibrium distribution functions 

fi^'") (e) = T ^ ia;„, - e " e-C^^-i)/^ exp(/3e) + 1 ' ^^^^ 
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A straightforward calculation of Z^^'*^) for the case of even A'' leads to the following expression 



^ ' iTTl{2k-l) 



5](-l)'-iexp(/?e(iV-0)exp(- 
i=\ ^ 



N 

^'"'w = SRAwn • 

where fc = 1, 7V/2. The equilibrium distribution functions (EDF) /(2S+i,fe) £qj. ^j^^g auxiliary Fermi-fields representing 
arbitrary S for SU{2) algebra are given by 

2S+1 



i-iy-' exp (/3e(25 + 1-0) exp 



-.(2S+i,fc)/ X ^ J^i ^ L (54) 

^ ^ ' exp((25 + l)/3e) + (-1)25+1 ^ 

for k = 1, [5 + 1/2J . Particularly simple are the cases of 5 = 1/2 and •5=1, 

= _ |„,(3,) _ iv^i^±f§ (55) 

2 2 smh(3pe/2j 

Here, the standard notations for Fermi/Bose distribution functions np/g{e) = [exp(/3e) ± 1]"-"^ are used. For 5=1/2 
the semi-fermionic EDF satisfies the obvious identity |/''^'^''(e)|2 = 71^(26). 

In general the EDF for half-integer and integer spins can be expressed in terms of Fermi and Bose EDF respectively. 
We note that since auxiliary Fermi fields introduced for the representation of SU(7V) generators do not represent the 
true quasiparticles of the problem, helping only to treat properly the constraint condition, the distribution functions 
for these objects in general do not have to be real functions. Nevertheless, one can prove that the imaginary part of 
the EDF does not affect the physical correlators and can be eliminated by introducing an infinitesimally small real 
part for the chemical potential. In spin problems, a uniform/staggered magnetic field usually plays the role of such 
real chemical potential for semi-fermions. 



II. APPLICATION OF SEMI-FERMIONIC REPRESENTATION 



In this section we illustrate some of the applications of SF representation for various problems of strongly correlated 
physics. 



A. Heisenberg model: FM, AFM and RVB 

The effective nonpolynomial action for Heisenberg model with ferromagnetic (FM) coupling has been investigated in 

[10]. The model with antiferromagnetic (AFM) interaction has been considered by means of scmi-fcrmionic represen- 
tation in [16] and [17] (magnon spectra) and in [11] for resonance valence bond (RVB) excitations. The Hamiltonian 
considered is given as 

Hint = - X! -^ii i^i^i ~\) (^^) 

<ij> ^ ' 



• Ferromagnetic coupling J = IpM > 

The exchange SiSj is represented as four-semi-fermion interaction. Applying the Hubbard-Stratonovich transforma- 
tion by the local vector field $ j (r) the effective nonpolynomial action is obtained in terms of vector c- field. The FM 
phase transition corresponds to the appearance at T < Tc of the nonzero average ($^(g = 0, 0)) which stands for the 
nonzero magnetization, or in other words, corresponds to the Bose condensation of the field 
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In one loop approximation the standard molecular field equation can be reproduced 

M = lFM{0)tanh{pM/2). 
The saddle point (mean-field) effective action is given by well-known expression 



Ao[M] = -N 



4/m(0) 



In 2 cosh 



I3M 



and the free energy per spin /o is determined by the standard equation: 

BM"^ ( 
/3/o = -InZs = i—- - In 2cosh 

Calculation of the second variation of >le// gives rise to the following expression 



(58) 
(59) 

(60) 



2cosh^(/3n) 



-^$+(fc,c.) 



tanh(/3Q) 



20 — iu) 



$-(A-,w) 



(61) 



where O = {gjisH + M)/2. The magnon spectrum (T < Tg) is determined by the poles of ) correlator, 

uj = Ak2. 

• Antiferromagnetic coupling J = Iafm < 0. Neel solution 
The AFM transition corresponds to formation of the staggered condensate 

$^(fc, w) = NmY'H^^^d^^o + ^'{k, oj) (62) 
The one-loop approximation leads to standard mean-field equations for the staggered magnetization 



= -Iafm{Q) tanh(/3Ar/2), Ao[N]=N 



4:Iafm{Q) 



In 2 cosh 



\ 2 J 



(63) 



After taking into account the second variation of Aef / , the following expression for the effective action is obtained 
[(see e.g. [16], [17]): 



^AFAlik) + 



+ ^$+(fc,a;) 



^AFAlik) + 



2cosh^(/3fi) 
2^2tanh(/3i7) 



$^(fc,0) + i lAFM{kWik,u^W{k,u^)- 



402 + ^2 



^- (fc, a;) - ^ $+(fc + Q, a;) $-(fc, a;). 



(64) 



The AFM magnon spectrum u = c|k|. 

• Antiferromagnetic coupling. Resonance Valence Bond solution 

The four-semi-fermion term in (56) is decoupled by bilocal scalar field Ajj . The RVB spin liquid (SL) instability in 
2D Heisenberg model corresponds to Bose-condensation of exciton-like [27] pairs of semi-fermions: 



Ao = - ^ ^ tanh (^) , Ao = - ^ In [2 cosh(/3/, Ao)] 

where Aq = A(q = 0) is determined by the modulus of A^- field 

A<ij>{R, f) = A(f)exp [ifA{R)^ 
whereas the second variation of SAej / describes the fluctuations of phase Aij 

Aeff=J2 ^"(k, u;)<f^A/3(k, c), <f^ = Tr{p"pf\Gj,+kGp + Gp+kGj,) + (^^^/(/pAo)) 

The spectrum of excitation in uniform SL is determined by zeros of tt^ and is purely diffusive [28]- [29]. 



(65) 

(66) 

(67) 
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B. Kondo lattices: competition between magnetic and Kondo correlations 



The problem of competition between Ruderman-Kittel-Kasuya-Yosida (RKKY) magnetic exchange and Kondo cor- 
relations is one of the most interesting problem of the heavy fermion physics. The recent experiments unambignonsly 
show, that such a competition is responsible for many unusual properties of the integer valent heavy fermion com- 
pounds e.g. quantum critical behavior, unusual antiferromagnetism and superconductivity (see references in [19]). 
We address the reader to the review [30] for details of complex physics of Kondo effect in heavy fermion compounds. 
In this section wc discuss the influence of Kondo effect on the competition between local (magnetic, spin glass) and 
non-local (RVB) correlations. The Ginzburg-Landau theory for nearly antiferromagnetic Kondo lattices has been 
constructed in [19] using the semi-fermion approach. We discuss the key results of this theory. 

The Hamiltonian of the Kondo lattice (KL) model is given by 

H = Y1 ^''^aCka + JJ2 (Sjsj + \Njnj) (68) 

fccr j ^ ^ 

Here the local electron and spin density operators for conduction electrons at site j are defined as 

Uj = J2 c]aCj'^ ' = ^ \(^]cTT'<ycr' Cja' , (69) 

j(T a 

where f arc the Pauli matrices and Cjo- = 5^^. Cka cxp{ikj). The spin glass (SG) freezing is possible if an additional 
quenched randomness of the inter-site exchange Iji between the localized spins arises. This disorder is described by 

H' = Y.Iji{SjSi). (70) 

We start with a perfect Kondo lattice. The spin correlations in KL are characterized by two energy scales, i.e., 
I ^ J'^/ep, and Ak ^ £f cxp(— ei?/J) (the inter-site indirect exchange of the RKKY type and the Kondo binding 
energy, respectively). At high enough temperature, the localized spins are weakly coupled with the electron Fermi 
sea having the Fermi energy e^, so that the magnetic response of a rare-earth sublattice of KL is of paramagnetic 
Curie- Weiss type. With decreasing temperature either a crossover to a strong-coupling Kondo singlet regime occurs 
at T ~ or the phase transition to an AFM state occurs at T = Tjv ~ zl where 2; is a coordination number in 
KL. If Tjv w Ax the interference between two trends results in the decrease of both characteristic temperatures or 
in suppressing one of them. The mechanism of suppression is based on the screening effect due to Kondo interaction. 
As we will show, the Kondo correlations screen the local order parameter, but leave nonlocal correlations intact. The 
mechanism of Kondo screening for single-impurity Kondo problem is illustrated on Fig. 5 




a) b) c) d) 

FIG. 5. Kondo screening of the local moment in single-impurity Kondo problem. Dashed line denotes semi-fermions, solid 
line stands for conduction electrons. 



As a result, the magnetization of local impurity in the presence of Kondo effect is determined in terms of GF's of 
semi-fermions G{lu) by the following expression [31]: 

M{H) = S{gnB)TY,{G^{'^) - Gli'^)) = S{gnB) ta.nh 

To take into account the screening effect in the lattice model we apply the semi-fermionic representation of spin 
operators. In accordance with the general path-integral approach to KL's, we first integrate over fast (electron) degrees 



1 - 



1 H\n{T/TK)) 
HT/Tk) 2\ti^{T/Tk) 



(71) 
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of freedom. The Kondo exchange interaction is decoupled by auxihary field (j) [32] with statistics complementary to 
that of semi-fermions which prevents this field from Bose condensation except at T = 0. As a result, we are left with 
an effective bosonic action describing low-energy properties of KL model at high T > Tk temperatures. 

• Kondo screening of the Neel order 

To analyze the influence of Kondo screening on formation of AFM order, we adopt the decoupling scheme for the 
Heiscnberg model discussed in Section II. A. Taking into account the classic part of Neel field, we calculate the 
Kondo-contribution to the effective action which depends on magnetic order parameter M: 



q,n 

where a polarization operator H{N) casts the form 

n(AA)=p(0)ln(^)- 



n(^) 



l<^n(q)l'- 



cosh(^AA) 



1 



-(C) 



(72) 



(73) 



where p{Q) is the density of states of conduction electrons at the Fermi level and the Kondo temperature Tk = 
eFexp(— l/(p(0)J)). Minimizing the effective action A{<I),M) with respect to classic field jV, the mean field equation 
for Neel transition is obtained (c.f. with (58)) 



N = tanh 



2T 



1 - 



atq coslr(J/Q.V72) 



In (T/Tk) cosh'(/37Q^) 



(74) 



As a result, Kondo corrections to the molecular field equation reduce the Neel temperature 

• Kondo enhancement of RVB correlations 

Applying the similar procedure to nonlocal RVB correlations, we take into account the influence of Kondo effect on 
RVB correlations 



n(7qA)=p(0)ln(^)+^ 



1 



cosh/?(/kA) 



- l + /kAtanh(/3/kA) 



k+q 



(75) 



Here = e(k) — ep- Minimizing the effective action with respect to A we obtain new self-consistent equation to 
determine the non-local semi-fermion correlator. 



tanh 



IqA 



/qA 



T\n{T/TK) 



(76) 



It is seen that unlike the case of local magnetic order, the Kondo scattering favors transition into the spin-liquid state, 
because the scattering means the involvement of the itinerant electron degrees of freedom into the spinon dynamics. 

• Kondo effect and quenched disorder 

Let's assume that the RKKY interactions are random (e.g. due to the presence of non-magnetic impurities resulting in 

appearance of random phase in the RKKY indirect exchange). In this case the spin glass phase should be considered. 
As it has been shown in [15] and [19], the influence of static disorder on Kondo effect in models with Ising exchange 
on fully connected lattices (Sherrington-Kirkpatrick model) can be taken into account by the mapping KL model 
with quenched disorder onto the single impurity Kondo model in random (depending on replicas) magnetic field. It 
allows for the self-consistent determination of the Edwards- Anderson g^^order parameter given by the following set 
of self-consistent equations 



2c 



\n{T/TK) 



O 



l + 2c{l3If{q-q)/\n{T/TK) 



O 



\n\T/TK] 



(77) 



Here q = qea and q are nondiagonal and diagonal elements of Parisi matrix respectively. Therefore, the Kondo- 
scattering results in the depression of the freezing temperature due to the screening effects in the same way as the 
magnetic moments and the one-site susceptibility are screened in the single-impurity Kondo problem (c.f. Fig. 5) when 
Ising and Kondo interactions are of the same order of magnitude. 
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FIG. 6. Feynman diagrams for nonlocal excitations associated with the overlap of Kondo clouds. 

Let's now briefly discuss the fluctuation efii'ects in Kondo lattices. The natural way to construct the fluctuation 
theory is to consider the non-local dynamical Kondo correlations described by the field 0(q, (sec Fig. 6). In fact, 
the non-locality of the "semi-Bosonic" field is associated with an overlap of Kondo clouds [19] and responsible for a 
crossover from the localized magnetism to the itinerant-like fluctuational spin-liquid magnetism. The temperature 
dependence of static magnetic susceptibility becomes nonuniversal in spite of the fact that we are in a region of critical 
AFM fluctuations which is consistent with recent experimental observations. 

C. Kondo efFect in quantum dots 

The single electron tunneling through the quantum dot [33] has been studied in great details during the recent 

decade. Among many interesting phenomena behind the unusual transport properties of mcsoscopic systems, the 
Kondo effect in quantum dots, recently observed experimentally, continues to attract an attention both of experimental 
and theoretical communities. The modern nanoscience technologies allow one to produce the highly controllable 
systems based on quantum dot devices and possessing many of properties of strongly correlated electron systems. 
The quantum dot in a semiconductor planar heterostructure is a confined few-electron system (see Fig. 7) contacted 
by sheets of two-dimensional gas (leads). Junctions between dot and leeds produce the exchange interaction between 
the spins of the dot and spins of itinerant 2D electron gas. Measuring the dc I — V characteristics, one can investigate 
the Kondo effect in quantum dots under various conditions. 

Various realizations of Kondo effect in quantum dots were proposed both theoretically and experimentally in recent 
publications (see e.g. [34] for review). In order to illustrate the application of semi-fermionic approach we discuss 
briefly electric fleld induced Kondo tunneling in double quantum dot (DQD). As was noticed in [35], quantum dots 
with even Af possess the dynamical symmetry 5*0(4) of spin rotator in the Kondo tunneling regime, provided the 
low-energy part of its spectrum is formed by a singlet-triplet (ST) pair, and all other excitations arc separated from the 
ST manifold by a gap noticeably exceeding the tunneling rate 7. A DQD with even A/" in a side-bound configuration 
where two wells are coupled by the tunneling v and only one of them (say, I) is coupled to metallic leads {L, R) is a 
simplest system satisfying this condition [35] . Such system was realized experimentally in Ref. [36] . 

Left lead (L) 
W 

Right lead (R) 

a) b) 
FIG. 7. (a) Double quantum dot in a side-bound configuration (b) co-tunneling processes in biased DQD responsible for the 
resonance Kondo tunneling. 

As it was shown in [20] the Shrieffer- Wolff (SW) transformation, when applied to a spin rotator results in the 
following effective spin Hamiltonian 

Hint = ^ Jaa'flfsc\o,„Ckia'cr+ ^ i^Jaa'^KA' + ^aa'Pth'^ '^ta'Aaa'^k''^''^' fhf^' C''^) 

kk' ,aa'—L^R kk',aa'AA' 
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where the c-operators describe the electrons in the leads and /-operators stand for the electrons in the dot. The 
matrices S"^ and P'^ {d~x,y,z) are 4x4 matrices defined by relations (43) (see Section I.C) and — J^^ , = J^^ 
and J^'^ are singlet, triplet and singlet-triplet coupling SW constants, respectively. 

Applying the semi-fermionic representation of S0(4) group introduced in Section I.C we started with perturbation 
theory results analyzing the most divergent Feynman diagrams for spin-rotator model [20] . Following the "poor man's 
scaling" approach we derive the system of coiiplcd renormalization group equations for effective couplings responsible 
for the transport through DQD. As a result, the differential conductance G{eV,T)/Go ~ \ Jlr\'^ shown to be the 
universal function of two parameters T/Tk and V/Tk, Go = c^/ttTi: 

G/Go ~ In"^ (max[(ey - 5), T]/Tk) (79) 

Thus, the tunneling through singlet DQDs with 5 = Et — Es ^ Tk exhibits a peak in differential conductance at 
eV K. 5 instead of the usual zero bias Kondo anomaly which arises in the opposite limit, 5 < Tk- Therefore, in 
this case the Kondo effect in DQD is induced by a strong external bias. The scaling equations can also be derived 
in Schwinger-Keldysh formalism (see [11] and also [18]) by applying the "poor man's scaling" approach directly to 
the dot conductance. The detailed analysis of the model (78) in a real-time formalism is a subject for a separate 
publication. 

III. EPILOGUE AND PERSPECTIVES 

In this paper, we demonstrated several examples of the applications of semi-fermionic representation to various 

problems of condensed matter physics. The list of these applications is not exhaustive. We did not discuss, e.g., the 
interesting development of SF approach for the Hubbard model with repulsive [14] and attractive [13] interaction, 
Dickc model, 2D Ising model in transverse magnetic field, application of SF formalism to mesoscopic physics [21] etc. 
Nevertheless, we would like to point out some problems of strongly correlated physics where the application of SF 
representation might be a promising alternative to existing field-theoretical methods. 

Heavy Fermions 

• Crossover from localized to itinerant magnetism in Kondo lattices 

• Quantum critical phenomena associated with competition between local and nonlocal correlations 

• Nonequilibrium spin liquids 

• Effects of spin impurities and defects in spin liquids 

• Crystalline Electric Field excitations in spin liquids 

• Dynamic theory of screening effects in Kondo spin glasses. 
Mesoscopic systems 

• Nonequilibrium Kondo effect in Quantum Dots 

• Two-channel Kondo in complex multiple dots 

• Spin chains, rings and ladders 

• Nonequilibrium spin transport in wires 

Summarizing, we constructed a general concept of semi-fermionic representation for SU(N) groups. The main ad- 
vantage of this representation in application to the strongly correlated systems in comparison with another methods 
is that the local constraint is taken into account exactly and the usual Feynman diagrammatic codex is applicable. 
The method proposed allows us to treat spins on the same footing as Fermi and Bose systems. The semi-fermionic 
approach can be helpful for the description of the quantum systems in the vicinity of a quantum phase transition 
point and for the nonequilibrium spin systems. 
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